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cn i ABSTRACT 

m . 

We revisit the feasibility of a cosmological test with the geometric distortion fo- 
' cusing on an ambiguous factor of the evolution of bias. Starting from defining es- 

\ timators for the spatial two-point correlation function and the power spectrum, in 

^ c"| a rigorous manner, we derive useful formulas for the two-point clustering statistics 

which take the light-cone effect and the redshift-space distortions into account. 

o ■ 

' Then we investigate how the predicted correlation functions are affected by the 

^ | redshift-space distortions and the bias assuming quasar samples which roughly 

^ , match the 2dF survey. Effect of the bias is crucial, in particular, on anisotropic 

•rH . 

■ component of the clustering statistics in redshift-space. By surveying behavior of 

5-h ' 

C$ ' the predicted correlation functions to parameters of a phenomenological model for 

the evolution of bias, it is shown that the correlation functions on the light-cone 
is sensitive to a mean amplitude of the bias and rather insensitive to the speed of 
its redshift-evolution. This feature is useful for an analysis of the geometric test. 
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1 INTRODUCTION 

The wide-field redshift survey projects, Two-degree Field (2dF) and Sloan Digital Sky Survey 
(SDSS), are upcoming, and large numbers of quasars and galaxies are expected to be detected. 
In particular, the 2dF group has recently reported their preliminary results about the evolution 
of the quasar luminosity function and the spatial correlation function based on the several 
thousands quasars detected from the survey so far (Shanks et al. 2000; Boyle et al. 2000). Q 
The SDSS survey project will also provide a set of high precision data for the distribution 
of quasars in near future. These surveys will provide very important clues on the origin and 
the evolution of the quasars, and will put constraints on theoretical models of the structure 
formation of the universe. However, there seems to remain room for discussion as to how we 
compare the observational data with theoretical models and what information is drawn out from 
them. Concerning the latter subject, various cosmological tests with the future data have been 
proposed. Among such cosmological tests, the geometric test using the cosmological redshift- 
space (geometric) distortion is quite unique. 

Alcock and Paczynski first pointed out the possibility of the geometric test (1979). They 
pointed out that if spherically distributed structures of cosmological objects were distributed 
in the universe, the observation of such structures at various redshift in redshift space offers a 
plausible test for the cosmological model of the universe, which is sensitive to a cosmological 
constant. This test is based on the fact that the shape of the spherically distributed structure 
is observed distorted in the cosmological redshift-space, where an incorrect distance-redshift 
relation is assumed. We refer to the distortion as the cosmological redshift-space (geometric) 
distortion. The geometric distortion is traced back to a coordinate transformation of the dis- 
tance between real space and redshift space. Several authors have extensively discussed about 
possible tests for cosmological models with the geometric distortion with the clustering statistics 
of high-redshift objects (Ryden 1995; Matsubara & Suto 1996; Ballinger et al. 1996; Popowski 
et al. 1998; Nair 1999; Hui et al. 1999). 

In principle, the geometric test is a unique test for the cosmological model, however, a 
measurement of the geometric distortion suffers from various other observational effects. In 
particular the peculiar motion of the cosmological objects causes additional distortion in the 

* After we have completed the present manuscript, the most recent result on the 2dF quasar survey has been reported (Croom et 
al. 2000), in which a clustering analysis of the quasars is presented in detail. 
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distribution in redshift space. Especially the linear distortion, which is due to the bulk motion 
of the cosmological objects, can be influential even at high-redshift on large scales. The net 
of the effect depends on the amplitude of the bias. Furthermore, the light-cone effect would 
become an important effect when the statistical quantities like the correlation function and 
the power spectrum are computed using the data in a wide range of the redshift. In this case, 
the evolution of the bias is an important but ambiguous factor in predicting the correlation 
function and the power spectrum (Suto, Magira, Yamamoto 2000; see also references therein). 
The theoretical studies based on the Press-Schechter theory and the extensive numerical works 
recently show that the time evolution of clustering bias for the galaxies and the clusters of 
galaxies strongly depends on the cosmological models (e.g, Mo & White 1996; Blanton et al. 
1999; Somerville et al. 2000; Taruya & Suto 2000), however, the quantitative prediction for the 
evolution of galaxy clustering is still under consideration. As for the quasars, the uncertainty 
of the formation mechanism makes for further difficulty in predicting the evolution of bias. 

In the present paper we consider the details of the theoretical predictions for quasar two- 
point statistics focusing on the evolution of the bias factor. We discuss how we can draw out 
information about the evolution of the bias from the two-point statistics of cosmological objects 
and how the evolution of bias alters the prediction for the two-point statistics, which determines 
feasibility of the cosmological test with the geometric distortion. This paper is organized as 
follows: In section 2, we summarize theoretical formulas for the two-point correlation function 
and the power spectrum which incorporate the light-cone effect and the redshift-space distortion 
effects. The final expression of the formula has been presented in reference (Suto, Magira, 
Yamamoto 2000). On the contrary to this previous paper, we here derive the theoretical formulas 
in a rigorous manner by defining estimators for the two-point correlation function and the power 
spectrum on a light cone. This makes clear how our theoretical two-point correlation function 
and power spectrum are related to the estimators in data processing. Then we investigate the 
details of the correlation functions assuming quasar samples which roughly match the 2dF 
survey in section 3. Section 4 is devoted to summary and discussions. Throughout this paper 
we use the unit in which the light velocity c equals 1. 



© 2000 RAS, MNRAS 000, 



4 



2 FORMALISM 

In this section we summarize formulas for the two-point statistics of the distribution of cos- 
mological objects in redshift space incorporating the light-cone effect. On the contrary to the 
previous paper by Suto, Magira, & Yamamoto (2000), in which the final expression for the two- 
point correlation function is presented, we here present a review to derive the formula starting 
from estimators for the anisotropic correlation function and the power spectrum in cosmological 
redshift-space. 

We start from reviewing an idealized data processing to compute the two-point correlation 
function and the power spectrum from a map of cosmological objects. Assuming a redshift 
survey, data is given in terms of a set of pairs of the direction vector 7 and the redshift z of the 
cosmological objects. When constructing the three dimensional map of the cosmological objects, 
we must introduce a radial coordinate s by assuming some distance-redshift relation s = s(z). 
Because the cosmological model of our universe is not definitely determined, we must assume 
the distance-redshift relation s = s(z). Frequently the distance extrapolating the Hubble law, 
s — s B.b( z ), or the comoving distance in the Einstein de Sitter universe, s = SEd( z ), have been 
used for the relation, where we defined 

SRb(z) = — , (1) 
-"0 



SEd(z) = — (1 - 1/VT+z), (2) 

with the Hubble constant Hq = 100/ikm/s/Mpc. We refer to the three dimensional space (5,7) 
as the cosmological redshift-space. For a sample at high-redshift z <S 1, extrapolation of the 
Hubble law, s = snb( z ), is unphysical. However, for comparison, we consider the both s = snb(z) 
and s = SEd(^) in sections 3.1 and 3.2, according to the previous investigations (Matsubara & 
Suto 1996; Ballinger et al. 1996). After section 3.3 we adopt s = SEd(^). 

Next we assume that the number density field n(s) with s = 57 can be constructed from 
the map. By introducing a synthetic catalogue n syn (s), we construct the density field: 

F(s) = -f-^) (3) 
[ / dsh(s) 2 ] 1 / 2 

where h(s) denotes the mean number density at the distance s. With the use of n(s), we write 
n(s) = n(s)(l + A(s)), (4) 



n 



syn 



s)=n( S )(l + A syn (s)). (5) 
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The synthetic field n syn (s) satisfies (ra syn (si)n syn (s 2 )) = (^(si)n syn (s 2 )) = n(si)n(s2), which 
lead 

(A syn ( Sl )A syn (s 2 )) = (A syn ( Sl )A(s 2 )) = 0. (6) 

Then the two-point correlation function can be computed by averaging the products of F(s) 
over all the possible pairs of two-points separated by the distance R. A mathematical expression 
for this process can be explicitly written 

&(R) = J ^ Jda 1 J ds 2 F(s 1 )F(s 2 )5^(s 1 - s 2 - R)^(s ft • R)(2Z + 1), (7) 

where R = |R|, R = H/R, s h = (si + s 2 )/2, s h = \s h \, s h = s h /sh, and Li(fi) is the Legendre 
polynomial of the Z-th order, (see figure 1.) In the case 1 = 0, where L (fj,) = 1, £q(R) represents 
the isotropic (angular averaged) correlation function. On the other hand, other cases measure 
the anisotropies of the correlation function owing to the factor Lj(s^ • R). In the particular case 
Z = 2, £ 2 (-R) expresses the quadru-pole anisotropy of the correlation function. We adopt this 
estimator for the isotropic and the anisotropic correlation functions. 

Similar to the case of the correlation function we define the following estimator for the 
power spectrum (cf., Yamamoto, Nishioka, & Suto 1999) 

Pi{k) = J ^ Jdsx J ds 2 F( Sl )F(s 2 )e ik -( Sl - S2 ^(s /t ■ k)(2Z + 1), (8) 

where k = |k|, k = k/Zc, and k is the wave number vector. When adopting the above estimators 
for the correlation function and the power spectrum, which we have defined independently, 
we can show that a familiar relation holds between them. Namely, by using the mathematical 
formulas, 

5 (3) ( Sl - s 2 - R) = J dke ik -( Sl - S2 - R \ (9) 

and 



e' kR 



A^ l ji(kR)Yr(^)Yr(n A ), (io) 

Lm 



where Y^ m (f2) and ji(x) denote the spherical harmonics on an unit sphere and the spherical 
Bessel function, respectively, we can verify that the following relation holds 



Zi(R) = 2^iiJ dkk^ikR^k). (11) 
The inverse transformation leads 

P t (k) = Am 1 J dRR 2 ji{kR)ii(R) . (12) 
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As is well recognized as the problem of the cosmic variance, expectation values of the 
estimators are only predictable theoretically. Then we consider the ensemble average of the 
estimators, i.e., 



MR)), 

(Pi(k)), 



(13) 
(14) 



which define the correlation function and the power spectrum on a light cone. 

In Appendix we show that, applying the small angle approximation, P^ c (k) is reduced to 

>dN\2/s 2 ds\ -i 

' h; '-'(A-..;) 

(15) 



/^(f) 2 (^)"V(M) 



r , /dN\ 2 ,s 2 ds\ -i 

J dz \~dV 

where dN/dz denotes the number count of the objects per unit redshift and per unit solid angle, 
Pr d (k,z) is defined: 

21 + 1 



C 2 Cm 



X 



d/j,Li(/j,)b(q[k], z) 2 P mass (q[k], z) 

i + / i 2 ((i + ^)v 2 -i; 



l + ^ 2 {uo 2 -l) 



2 D \ jika p (z) 



(16) 



where P mass (q, z) denotes the power spectrum of the total mass distribution in real space with 
the wave number q, D[fika p (z)/c n ] describes the damping function due to the Finger of God 
effect, and with denoting the comoving radial distance in real space r(z) and the distance in 
cosmological redshift-space s(z) (see equations (1), (2) and (A27) for the definitions), we defined 

r(z) 



S ( z y 

dr(z) 
ds(z) 



q[k] = k 



\ 



;i-m 2 ) , ^ 

c 2 c 2 , ' 



and 



i 



dlnD^z) 



(17) 

(18) 
(19) 

(20) 
(21) 



b(q[k],z) d\na(z) 

with the linear growth rate Di(z) and the scale factor a(z). The function b(q,z) denotes the 
linear bias factor defined in real space (see Appendix). Expression (|TTD is correct only for the case 
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when the real space is spatially flat, (see equation ( |A28|) for generalization to an open universe) 
The correlation function ^ LC (i?) can be computed from Pj jC {k) with the use of equation (|TT|). 

The expression P z crd (/c, z) is well-known as a power spectrum on a constant time hypersurface 
in the cosmological redshift-space (e.g., Ballinger et al. 1996; Suto et al. 1999; Magira, Jing, 
& Suto 2000). The term including the /3(z)-factor in the expression P[ rd (k,z) is originated 
from the Kaiser factor which represents the effect of the linear distortion (Kaiser 1987). In 
the present paper we adopt the exponential model for the distribution of the pairwise peculiar 
velocity, then the damping factor is 

D| " to ' (2)/C ' 1 = rr^wv^ ' (22) 

where a p (z) is the pairwise velocity dispersion. For P maaa (q,z), the power spectrum in real 
space at the redshift z, we adopt the linear power spectrum with the CDM transfer function 
(Bardeen et al. 1986) and the Peacock & Dodds formula for the nonlinear power spectrum 
(Peacock & Dodds 1994; 1996). The essence of the geometric distortion is only the coordinate 
transformation between (s(z), 7) and (r(^),7), which is described by the coefficients c ± and Cy. 

Equation ( [Op indicates that the light-cone effect is incorporated by averaging the local 
power spectrum P" d (k, z) by weighting the factor in the formula: 

*w-(£)v=r- <»> 



The formula fllq) is obtained under the small angle approximation. Hence the validity of the 
use is limited to small scales. For large scales, R *S s max , where s max is the size of a survey area, 
the finite size effect of the survey area affects the proper estimation of the correlation function 
and the power spectrum, (e.g., Yamamoto & Suto 1999; Nishioka & Yamamoto 1999; 2000). 



3 QSO TWO-POINT STATISTICS AND THE EVOLUTION OF BIAS 

In this section we discuss the detail of the two-point correlation function focusing on how the 
quasar two-point statistics depend on the evolution of bias. Motivated from the recent report 
by Shanks et al. (2000), we assume a sample which roughly corresponds with the 2dF quasar 
survey. 
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3.1 The Quasar Luminosity Function 

For simplicity we here adopt the B-band quasar luminosity function according to Wallington & 
Narayan (1993). We set the limiting magnitude B\ [m = 20.85, which is useful to obtain a sample 
which roughly matches the 2dF quasar survey. According to Wallington & Narayan (1993; see 
also Nakamura & Suto 1997), we adopt 

$(M B ,z) = $* X [10°- 4 [ M B-iW B (2)](«i+l) _|_ 1 Q0.4[M B -Afe(*)]('e2+l)]-l j (24) 

where 

$* = 6.4 x 10~ 6 h 3 Mpc- 3 mag _1 , (25) 

M B (z)=M*-2.5k L log(l + z), (26) 

M B = —20.91 + 5 log h, (27) 

with ki = 3.15, K\ = —3.79, = —1.44, which is applicable for 0.3 < z < 2.2. The B-band 
apparent magnitude is computed from a quasar of absolute magnitude M B at z, 
~d L (z) 



B = M B + 5 log 



-2.5(1 -p)log(l +z), (2* 



lOpc 

where cIl(z) is the luminosity distance and we adopt the quasar energy spectrum L(y) oc 
v~ v (p=0.5) for the K-correction. Because the luminosity function is constructed under the 
assumption of the Einstein de Sitter universe, we here use the luminosity distance 



d L (z) = — (1 + z- VT+z), (29) 

which holds in the Einstein de Sitter universe. 

The number count of the sample with the B-band limiting magnitude Bn m per unit redshift 
and unit solid angle is expressed 

dM B $(M B ,z), (30) 

az az Jm b ,. 

with r being the comoving distance. Because of the same reason for the luminosity distance 
d,L, we adopt the comoving distance in the Einstein de Sitter universe being r = cIl(z)/(1 + z). 
Figure 2 plots the weight factor W(z) as the function of redshift with adopting cosmo logical 
redshift-space extrapolating the Hubble law s = snb(z) (dotted line) and Einstein de Sitter 
universe s = SEd{z) (solid line). From this figure we see that the choice of the cosmological 
redshift-space significantly alters the weight factor. After completing almost of the present 
paper, the preliminary report on the luminosity function from the 2dF quasar survey (Boyle 
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et al. 2000) has been announced. They have found the consistency of the luminosity function 
with the previous results. The use of the new luminosity function does not alter our conclusions 
described below. 

3.2 QSO Spatial Clustering and Evolution of Bias 

Recently bias models for the quasar distribution have been proposed (e.g., Martini & Weinberg 
2000; Haiman & Hui 2000; Fang & Jing 1998). Construction of theoretical models for the quasar 
bias is a challenging problem, and we need more investigations. However, the bias mechanism 
of the quasar distribution contains many uncertain factors at present. In the present paper we 
consider one of the phenomenologically simplest models of the bias with the form, 

b(z) = a + {b{z m )- a )(j±^f, (31) 

where a, /3, and b(z*) are the constant free-parameters (cf., Matarrese et al. 1997). Throughout 
the present paper we use f3 to denote the free parameter of the bias evolution, which should 
be distinguished from (3(z) which we use to denote /3(z)-factor for the linear distortion. In 
equation (pi]), (3 specifies the rate of evolution of bias, b(z*) does the amplitude of bias at some 
redshift z*, and a corresponds to the amplitude of bias at z = in the limit 2* ^> 1. We 
have assumed the bias model is independent of scales k. This assumption relies on the recent 
numerical simulations that the scale-independent bias can be valid on the linear and quasi- 
linear scales as long as the mechanism of bias is dominantly affected by the gravity (e.g, Mann, 
Peacock & Heavens 1998). Although the realistic bias model taking account of the nonlinearity 
and stochasticity should be considered correctly (Dekel & Lahav 1999; Fry 1996; Tegmark & 
Peebles 1998; Taruya, Koyama & Soda 1999; Taruya & Suto 2000), the deterministic linear 
bias parameterized by the one parameter b can become a good approximation even in the 
quasi-linear regime for galaxies and clusters (Taruya et al. 2001; see also Scherrer & Weinberg 
1998). 

As noted above, both of the power spectrum Pf jC {k) and the correlation function ^ LC (i?), in 
principle, carry the same information because they are connected by the Fourier transformation. 
As the light-cone effect can be regarded as an averaging process, therefore, it is important to 
understand the evolution of P l crd (k,z) (or £f rd (R, z)) at different redshift. 

Figure 3 shows typical behaviors of the power spectrum PQ rd (A;, z) and P^ik, z) at dif- 
ferent redshifts adopting the cosmological redshift-space with s = SHb(z). The upper pan- 
es) 2000 RAS, MNRAS 000, [l|-§l] 
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els show PQ rd (k,z), the middle panels show P| rd (/c,z), and the lower panels show the ratio 
P| rd (/c, z)/PQ rd (k, z). The left panels show the case for z = 0.5, the center for z = 2.0, and the 
right for z — 4. Here we adopted the bias model given by equation ([H]) with a = 0.5, /3 = 
1, b(z* = 2) = 2, and the real space of the CDM model with a cosmological constant, to be spe- 
cific, Q = 0.3, Q\ = 0.7, h = 0.7, cr s = 1.0. The ACDM model reproduces the observed cluster 
abundance (Kitayama &; Suto 1997). Figure 4 is same as figure 3, but is the case adopting the 
cosmological redshift-space with s = SEd{z). 

In each panel, the dashed line plots the case when the distortion due to the peculiar motion 
of the sources are neglected, to be specific, the Kaiser factor, the nonlinear and the finger of 
god effect are neglected in the formula. Hence the dashed line shows the case when only the 
geometric distortion is taken into account. The dotted line shows the case when the nonlinear 
and the finger of God effects are neglected, and the solid line plots the case obtained from 
the formula (15), in which all the effects are incorporated. These two figures 3 and 4 indicate 
how the three kinds of the distortion effects, i.e., the geometric distortion, the linear distortion, 
the finger of god effect, are effective at different redshifts and on different scales. The three 
lines, the dashed line, the dotted line, and the solid line, in upper panels behave in a similar 
way. Therefore, the effect from the peculiar motion is small for the amplitude Pg rd (/c, z). The 
geometric distortion significantly affects the amplitude of the power spectrum Pg rd (fc, z). Note 
also that the effect acts in different ways for the cases s = SHb{z) and s = SEd{ z )- As the 
redshift becomes higher, the amplitude of Pg rd (fc, z) becomes larger (smaller) in the case with 
s = sm (s = SEd)- This difference is caused by the scaling effect through the coefficients c ± and 
C|| (See also Taruya & Yamamoto 2000). Thus, for the power spectrum, PQ rd (fc, z), the geometric 
distortion and the bias are the most important effects for the amplitude. 

We read the following behaviors from the middle and lower panels of figures 3 and 4. For the 
anisotropic power spectrum, the redshift-dependence is remarkable when adopting s = snb{ z )- 
In this case, the geometric distortion is a minor effect, however, the linear distortion and the 
finger of god effect are effective on the large scales and on the small scales, respectively, at 
low redshifts. At the higher redshifts, the geometric distortion also becomes influential. On the 
contrary to this, in the case adopting s = SEd(z), the geometric distortion effect is a minor 
effect compared with the linear distortion for P| rd (/c, z) even at high-redshift z = 4. Since we 
here fixed the bias parameter, this feature may be strongly dependent on the behavior of the 
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bias parameters. Therefore we investigate the bias-dependence of the correlation function in 
detail in the next subsections. 

From figures 3 and 4 it is apparent that the choice of the cosmological redshift-space affects 
the behavior of the power spectrum (the correlation functions) . To avoid the unphysical increase 
of the clustering amplitude at high-redshift in the case s = snb{z) and to compare with the 
observational result (Shanks et al. 2000), in which QSO correlation function is reported by 
assuming the Einstein de Sitter universe, hereafter we adopt s = SEd{z)- Though the behavior 
of the correlation functions depends on the choice of the cosmological redshift-space, sensitivity 
of the ratio of the correlation functions on a light-cone to the cosmological constant will not be 
significantly altered (see subsections 3.3 and 3.4), neither will be the feasibility of the geometric 
test. 

3.3 Characteristic Correlation Length Determines QSO Bias 

In this subsection we consider how the amplitude of the correlation functions which incorporate 
the light-cone effect depends on the evolution of bias. Here we determined of the bias model 
by the formula 

_ JdzzW(z) 

Z *~ JdzW(z) ' {6Z) 
which describes a mean redshift. Assuming the quasar sample with Bn m = 20.85 in the range 

0.3 < z < 2.3, we have = 1.3 adopting s = SEd(z)- We adopt this set of the values below. 

Now we define the characteristic correlation length R c by 

£ LC (i? c ) = 1. (33) 

Figure 5 shows contour of R c on the plane of the bias parameters b(z*) and /3, which is computed 
from the definition ( |3~3"D with the theoretical correlation function £q C . Hereafter, unless otherwise 
stated, the correlation functions are computed by taking into account all the effects, i.e., the 
linear distortion, the geometric distortion, the nonlinear and the light-cone effects described 
in section 2. The left panels show the case a = 0, the center a = 0.5, and the right a = 1.0, 
respectively. In the upper panels, the cosmological model with a cosmological constant, Qq — 
0.3, Qa = 0.7, A = 0.7, <7 8 = 1 is assumed as the model for real space, and the Einstein de 
Sitter universe is adopted as the cosmological redshift-space. The lower panels show the case 
assuming the real space of an open universe Q = 0.3, Qk = 0.7, h = 0.7, cr 8 = 1. It is notable 
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that R c does not significantly depend on the parameters a and ft, but is sensitive to b(z*), the 
amplitude of the bias at z*. This feature is caused because z* is roughly the mean value of the 
redshift for averaging the correlation function. 

Shanks et al. (2000) have reported that the QSO correlation function is consistent with being 
(r/r ) -1 ' 8 with r = 4/i _1 Mpc. This results is consistent with the previous paper (Croom & 
Shanks 1996), in which the correlation function was fitted in the similar form with same power 
index with ro = 6/i~ 1 Mpc. At present it seems that there exist not small statistical errors in 
the observed correlation function of the quasar distribution. If the value of R c was fixed, it will 
put a significant constraint on the evolution of bias. Though we need to assume the value of 
erg, the characteristic correlation length determines the QSO bias. For example, if we seriously 
take the result r = 4 ~ 6/i _1 Mpc by Shanks et al. (2000) and Croom & Shanks (1996), figure 
5 suggests that the bias of the quasar sample has the amplitude of factor 1 ~ 2 at = 1.3. 
However, this statement might not be taken so seriously, because the result is preliminary and 
the sample is not homogeneous at present. And also note that the estimators to evaluate the 
correlation function should be taken in a consistent way in data processing according to the 
theoretical computations. 

3.4 Feasibility of the Geometric Test 

The geometric distortion will be detected by measuring ^(-R), so we next consider the anisotropic 
part of the correlation function, C,2 C (R). The anisotropic power spectrum P 2 LC (A;) can also mea- 
sure the geometric distortion (Ballinger et al. 1996). In the present paper, however, we consider 
the ratio & C {R)/$ C (R) (cf., Matsubara & Suto 1996). Fi srure 6 shows the contour of the ratio 
£,2° (R) / (,q C (R) on the plane of the bias parameters b(z*) and ft, with fixed R = 5/i _1 Mpc 
(upper panels), R = 10/i _1 Mpc (middle panels), and R = 20/i _1 Mpc (lower panels). The left 
panels show the case a = 0, the center a = 0.5, and the right a = 1.0. In this figure the 
ACDM model is assumed as the model of real space, where the same cosmological parameters 
are taken as those in figure 5, and the Einstein de Sitter universe is chosen for the cosmological 
redshift-space. It can be read from this figure that the dependence on a is weak, and that the 
dependence on ft becomes rather weak for R ~ 10/i _1 Mpc. The amplitude of ^2 C /^o C ^ s m °st 
sensitive to the parameter b(z*). The reason would be same as that described in subsection 3.3 
for contour of R c . 
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To show the importance of the linear distortion, we show the contour of the ratio £,2 C /£o° 
with fixing R = 10/i _1 Mpc, in figure 7. In this figure the same model parameters as those in 
figure 6 are adopted, however, the distortion effects due to the peculiar motions are omitted on 
purpose in the upper and the middle panels. The upper panels show the case only the geometric 
distortion is incorporated but the distortion effect due to the peculiar motion of sources, i.e., the 
linear distortion, nonlinear and the finger of god effects, are neglected. The middle panels show 
the prediction within linear theory, that is, the case the linear distortion effect is incorporated 
but the nonlinear and the finger of God effects are omitted in the computation. The lower 
panels show the case all the effects are incorporated (Same as the middle panels in figure 6). 
From figure 7 it is apparent that the middle and the lower panels are significantly different from 
the upper panels. This means that the linear distortion, the nonlinear and the finger of god 
effects are the dominant effects for £g C a t this scale. In a redshift survey, the peculiar motion of 
sources inevitably contaminate the map. These facts limit the ability of the cosmological test 
with the geometric distortion because the linear distortion effect is sensitive to the amplitude 
of bias. 

Nevertheless we would like to discuss that the cosmological test with the geometric distortion 
might be a useful tool for testing the cosmological model of our universe with the precise 
measurement of the correlation functions. In order to discuss the feasibility of the geometric 
test, we also consider the case that the real space is an open hyperbolic universe. Prescription 
to compute the correlation functions for the open universe is described in Appendix. Figure 8 
shows the contour of £2° (R) / £,q C (R) to demonstrate the feasibility of the geometric test. The 
solid lines show contour of £2 (^)/£o (^-) with fixed R = 10/i _1 Mpc for various cosmological 
models (a) ACDM model with Q = 0.3; (b) ACDM model with Q = 0.4; (c) open CDM 
model with Q = 0.3; (d) open CDM model with Q = 0.4. In each model we adopted h = 0.7, 
a = 0.5, and as normalized by cluster abundance (Kitayama &; Suto 1997). The region between 
the dashed lines in each panel satisfies 4 < R c < 6 for each model. Therefore ^2 C /^o° a ^ 
R = 10/i~ 1 Mpc is in the range —0.8 ~ —0.5 for the ACDM model taking the range of the bias 
parameter consistent with R c = 4 ~ 6/i -1 Mpc into account. While ^2 C /^o° ~ —0.5 ~ —0.3 is 
predicted in the open CDM models. In this figure we fixed a = 0.5, however, the dependence 
on a is weak for ~ a ~ 1. 

Figure 9 shows the same figure 8 but with fixed R = 20/i~ 1 Mpc. In this case we find that 
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C2 C /Co° is in the ran g e -2-0 1.2 (-2.0 1.4) for the ACDM with Q = 0.3 (O = 0.4) 

taking the range of the bias parameter consistent with R c = 4 ~ 6^i _1 Mpc into account. While 

^2 C /^o° ~ -1-2 0.8 (-1.4 ~ -1.0) is predicted in the open CDM model with Q = 0.3 

(f^o = 0.4). This figure demonstrates that the difference in ^2 C /^o C between the ACDM model 
and the open model appears irrespective of the scale R. At the same time this figure shows that 
the predictions slightly depend on the value of Qq. If the constraint on R c becomes tight, the 
difference in ^2°/^o C becomes remarkable. The ratio £2 C /£o° does not depend on cr 8 at large 
length scales unless the nonlinear effect becomes influential. However, we should remind that the 
assumption of the value of as is needed to constrain the bias parameters by the characteristic 
correlation length. 

This difference in ^2 C /^o C between the ACDM model and the open model occurs due to the 
difference of the linear growth rate and the geometric distortion effect. The scaling effect due 
to the geometric distortion alters the amplitude of £q C (-R) through the coefficient c ± and Cy. 
The evolution of c± and Cy in the ACDM model is quite different from those in the open model. 
The main reason for the difference in ^2 C /^o C ^ s fr° m the difference in £q C due to the scaling 
effect of the geometric distortion. The difference appears for wide range of the parameters of 
the bias evolution. In order to check the model dependence of the bias, we have computed the 
ratio £,2 C /£o C with the following parameterization for the bias instead of (|3~ID , 



This alternation does not change the result significantly. These situation suggest that precise 
measurements of correlation functions might offer a unique and useful test of the cosmolog- 
ical model of our universe, however, further investigation should be needed for the definite 
conclusion. 



In the present paper we have revisited the feasibility of the geometric distortion focusing on 
the ambiguous factor of the evolution of bias. We have derived the useful formulas for the 
correlation function and the power spectrum which incorporate the light-cone effect and the 
redshift-space distortions in the rigorous manner by defining the estimators for the anisotropic 
correlation function and the power spectrum. The final expressions are consistent with those 




(34) 



4 SUMMARY AND CONCLUSIONS 
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presented in the previous paper by Suto, Magira, and Yamamoto (2000). Our investigation 
makes clear how the theoretical two-point correlation function and the power spectrum in our 
paper correspond to estimators in data processing. 

For a quasar sample like the 2dF survey, whose distribution is extended to high-redshift, the 
light-cone effect becomes very important. In the correlation functions, the light-cone effect is 
incorporated by averaging a local correlation function over the redshift with the weight factor 
W(z). With this formula we have investigated the theoretical predictions assuming a simple 
model of quasar sample which roughly match with the 2dF survey. Adopting a model for the 
evolution of bias ([H]), which is phenomenologically parameterized in the simple form, we have 
examined the theoretical predictions for the two-point correlation functions. As pointed out by 
Ballinger et al.(1996), the linear distortion is the dominant effect even at the high-redshift on 
the anisotropic correlation function £2°. The linear distortion is sensitive to the amplitude of 
bias, which limits the feasibility of the cosmological test with the geometric distortion. 

We have shown that the predicted correlation functions on the light-cone are sensitive to the 
amplitude of bias b(z*) at the mean redshift z*, and rather insensitive to the other parameters 
which specify the speed of its redshift-e volution due to the averaging process from the light-cone 
effect. The characteristic correlation length of the isotropic correlation function can determine 
the mean amplitude of the bias, though we need to assume the value of erg. If we take the 
discussions in subsection 3.3 seriously, the quasar correlation function is consistent with the 
ACDM model with the amplitude of the bias being 1.0 ~ 2.0 at z* = 1.3. We have also 
shown that ^ c (R) / C,q C (R) is only sensitive to the mean amplitude of bias, not to the speed 
of its redshift-evolution. Taking the constraint on the mean amplitude of the bias from the 
characteristic correlation length into account, the measurement of £,2(R)/£,o(R) can be a useful 
tool for testing the cosmological model of the universe. Conversely if the cosmological model is 
fixed, the precise measurement of £i(R) and £o(R) will provide a clue for the bias model. 

Finally we mention about other possible effects which affect measurements of statistical 
quantities. In general, the sampling noise and the shot noise limit the precise measurements 
of the statistical quantities. The shape of a survey area might affect a proper estimation of 
the correlation functions too. Those effects depend on observational strategies. We cannot 
conclude whether the 2dF quasar sample will detect the cosmological constant at present stage 
and a numerical approach is needed to examine those observational effects to draw definite 
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conclusions. Furthermore, though we have assumed the time-independent cosmological constant 
and the deterministic bias, the possibilities of the time-dependent cosmological constant like 
the quintessential CDM model and the stochasticity of the QSO clustering bias would make the 
feasibility of the geometric test more complex (Yamamoto & Nishioka 2000). N-body numerical 
simulations of a cosmological horizon size, which is the frontier of numerical cosmology (e.g., 
Pearce et al. 2000; Hamana et al. 2000), is another useful approach to investigate the geometric 
test. 
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APPENDIX A: 

In this Appendix we give an explicit derivation of the formula (15). Though the formula is 



presented in the paper (Suto et al. 2000), the explicit derivation is not presented there. We 
believe that it is useful to present the way to derive the formula from the estimators Pi(k) and 
£z(-R) defined explicitly. From the definition of the estimators, we write 

$ C (R) = Jds 1 Jds 2 (F( Sl )F( S2 ))8^\s 1 -s 2 -R)L l (s h -R)(2l + l). (Al) 



With the aid of equation (F(si)F(s 2 )) in ( |A1| ) reduces to 



(F(Sl)F(S2)) = Jd^sf • (A2) 

Instead of the variables s x and s 2 , we introduce 

x = Si - s 2 , (A3) 
s h = ^( Sl + s 2 ), (A4) 

then we have 

$ C (R) = [4tt J dss^s) 2 ]- 1 1 ^j d *j ds h n( Sl )n(s 2 )(A( Sl )A(s 2 )) 

x 5 (3) (x-R)L ; (s ft -R)(2/ + l), (A5) 

where Sh = |s/j| and Sh = s h /\s h \. Exact calculations for £q C (R) and P LC (/c) are performed 
within linear theory of density perturbations in references (Suto & Yamamoto 1999; Nishioka & 
Yamamoto 1999; 2000). Though all the redshift distortions are not incorporated simultaneously 
in those works, however, those investigations have shown the validity for using the distant 
observer approximation. Then we here use the distant-observer approximation: 

w(si)n(s 2 ) ~ n(s h ) 2 } (A6) 
(A(si)A(s 2 ))^^(s&,x,s /l -x), (A7) 

where x = x/|x| and x = |x|. Here we have approximated that (A(s;l)A(s 2 )) is a function of 
Sh, x, and • x. 
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Now we consider the expression for £(sh,x,Sh • x). We introduce the following coordinates 
to describe real universe, 



ds 2 = a^f^-drf + dr 2 + f K (r) 2 dn 2 2) ) , 
where 

' s\n{VKr)/VK (K > 0) 

f K (r) = \ r {K = 0) ; 



(A? 



(A9) 



H%Qk- Then the number conservation 



^ sinh (^Kr)fV^K (K < 0) 
with the spatial curvature K = Hq(Q + — 1) = 
means 

n(s, -f)s 2 ds = n R (r, ~/)f K (r) 2 dr, (A10) 

where n R (r, 7) denotes the number density field in real space. Introducing the mean number 
density n(r) R at the distance r in real space, we write 

n R (r,7) = n R (r)(l + A R (r,7)), (All) 

where A R (r, 7) denotes the density contrast in real space. Combining (|A10 ) and (|A11 ), we have 

n(s)s 2 ds = n K (r)f K (r) 2 dr, (A12) 

n(s) A(s, -f)s 2 ds = n R (r) A R (r, j)f K (r) 2 dr, (A13) 

and 

(A( Sl )A(s 2 )) = (A R ( ri )A R (r 2 )), (A14) 

where r denotes the point specified by r and 7. When the distance |r x — r 2 | is small compared 
with the size of the survey volume, we may write 

(A R ( ri )A R (r 2 )) 



dq^j 
(2tt) 3 



exp[zgu(ri - r 2 ) + iq ± ■ (r ±)1 - r ±>2 )]P{q {l , q x , z) 



(A15) 



where P(q h q ± , z) denotes the power spectrum of cosmological objects in real space (including 
the effect of the peculiar motion) on a constant time hypersurface at redshift z, and we intro- 
duced r = (r, r ± ), q = (q n , q ± ), where the subscript 1 1 denotes the component of the line of sight 
and _L denotes the components perpendicular to the vector of the line of sight. Here the redshift 
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z in P(q\\, q±, z) should be understood as a function of the comoving distance = (n + ^)/2. 
Even in the case of the open universe, we can expand the correlation function in terms of the 
usual Fourier modes, as in (|A15 ), because we are considering the case the distance |ri — r2 1 
is sufficiently small compared with the horizon (curvature) scale. Under the distant observer 
approximation we have 

r_L,i - r ± , 2 = lz^{a ±A - s ± , 2 ), (A16) 

Sh 

n -r 2 = — (si - s 2 ), (A17) 
where = |si + S2I/2. Introducing the new variables, 

, fK\Th) ( Kio\ 

k ± = q ± = c x q ± , (A18) 

Sh 

fcli = -^r<l\\ = C \\Q\\, ( A19 ) 



we have 



= ik-(si- S2 ) n/ ^11 k ± 



(A( Sl )A(s 2 )> = / e *^-«ip( g _> _JL, q± _> (A20) 
c^Cy 7 (27r) ,:i C|| c ± 

where we used ( |A15| ) and ( |A14|) . Here the redshift z should be understood as a function of Sh 
instead of r^. Inserting this equation into (|7A5|), we have 

£ C (R) = [4tt I daM'fT^J dn Sh jds h sln(s h ) 2 ^ J j^y 3 Ji(kR) 

xP(q n - h q± _> ^£ z ) Ll (g ft ■ k)(2/ + 1), (A21) 
where we used the mathematical formulas, 

e^ = 47r^ ]T ^(^)yr(^k)^ m (^x)> (A22) 

l m=—l 

Li(s h ■ k) = — — £ rr(^J^ m *(^ k ). (A23) 

m=—l 

By introducing /x = we have 

£ C (R) = [Jdss*h{s)r l i l J ds h sln{s h f 2 -^ f^-JiikR) £dp 

xP(g„ -> ^,|q ± | - AVr ^ ,z)L ; (/i). (A24) 

C11 c x 



The transformation (|12|) leads 

P Z LC (A;) = [[dss 2 h(s)T 1 I ds h s 2 h n{s h f 2 -^- C d^L^) 
J J c,c» Jo 
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xP( g|| ^,|q ± H^L^!,). (A25) 

Gii Cj_ 

The effects of the linear distortion, the nonlinear effect, and the finger of god effect have 
been investigated for the two-point statistics on a constant time hypersurface (Peacock & Dodds 
1994; 1996, Ballinger et al. 1996, Magira, Jing, & Suto 2000). We adopt the following formula, 



P(q h \q ± \,z) = {l + P(z)^)yb(q,z) 2 P mass (q,z)D[q l]C r P (z)], (A26) 



with q = Jq 2 + |q±| 2 , mass power spectrum P mass (q, z), the bias factor b(q,z), the damping 



factor D[q }l cr P (z)}, and the (3(z) factor defined by equation (|2l|) . Here we have assumed that 

the number density fluctuation of cosmological objects is simply proportional to the mass 

density fluctuation in real space by introducing the linear bias factor b(q,z). Then we obtain 

the expression (|i~5l) , using the relation dN = dss 2 n(s). 

For convenience we summarize the explicit formulas for the factors Cy and c x . Throughout 

this paper we assume that the correct model of our universe is the CDM model. Then the 

comoving distance in the real space is expressed as 

, s 1 f z dz' 
r{z) = — / A27 

^0^0 ^(1 + ^)3 + ^(1+^)2 + ^ 



With this formula, c ± is given by 

c ± = ^++\ < A28 > 
s(z) 

and 

' (A29) 



, (A30) 

^o(1 + ^) 3 + ^(1 + ^) 2 + ^a 
for the case s = SHb(z), respectively. 

This paper has been produced using the Royal Astronomical Society/Blackwell Science IsTgX 
style file. 



'n + n K (i + z)~ l + n A (i + z )- 3 ' 

for the case s = SEd( z ), and 
1 
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FIGURE CAPTIONS 

Fig. 1 — A sketch to illustrate variables defined in this paper. 

Fig. 2 — The weight function W(z) for the quasar sample in section 3.1. The dashed line shows 
the case adopting s = snb(z), and the solid line does the case s = SEd( z )- 

Fig. 3 — Power spectra on a constant-time hypersurface. The upper panels show PQ rd (k,z) 
in unit of /?r 3 Mpc 3 , the middle panels show P 2 crd (fc, z), and the lower panels show the ratio 
P 2 crd (fc, z)/PQ rd (k, z). The left panels show the case z = 0.5, the center z — 2, and the right 
z = 4. We assumed the ACDM model with fi = 0.3, fi A = 0.7, h = 0.7, and a s = 1.0. In this 
figure the bias parameter is fixed as a = 0.5, (3 = 1, and b(z* = 2) = 2. Here s = SHb{z) is 
adopted. In each panel, the solid line plots the case obtained from the formula (15), while the 
dotted line plots the case when the nonlinear and the finger of God effects are neglected, and 
the dashed line plots the case when the distortion due to the peculiar motion of the sources 
and the nonlinear effect are neglected. The dashed line shows the case that only the geometric 
distortion is taken into account. 

Fig. 4 — Same as Fig. 3 but with adopting the cosmological redshift-space s = SEd(z). 

Fig. 5 — Contour of the characteristic correlation length R c on the bias-parameter plane f3 and 
b(z*) with z* = 1.3. We fixed a = in the left panels, a = 0.5 in the center panels, and 
a = 1.0 in the right panels. The upper panels show the case when the ACDM model, Qq — 0.3, 
Q\ = 0.7, h = 0.7, cr 8 = 1.0, is assumed as the cosmological model of the real space. While 
the lower panels show the case that the open CDM model with Qq = 0.3, = 0.0, h = 0.7, 
ct 8 = 1.0, is adopted as the real space. We adopted the cosmological redshift space s = se<i{z)- 
In each panel, the lines show the contour of the levels, R c = 10, 8, 6, 4, 3 /i _1 Mpc, from right 
to left. The solid lines are the contours R c = 6 and 4 /i _1 Mpc. 

Fig. 6— Contour of & C (R)/^ C (R) in the ACDM model with fixed R = 5/r 1 Mpc (upper 
panels), R = 10/i~ 1 Mpc (middle panels), and R = 20/i~ 1 Mpc (lower panels). Similar to figure 
5, we chose a = (left panels), 0.5 (center panels), and 1.0 (right panels). In this figure all 
the distortion effects are taken into account. The cosmological parameters, f2 = 0.3, Q\ = 0.7, 
h = 0.7, as = 1.0, are taken, and we adopted the cosmological redshift space s = SEd{z). 
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Fig. 7 — Contour of £2 C (-^)/£o°(-^) with fixed R = 10/i _1 Mpc on the bias-parameter plane 
P and 6(2*). We fixed a = (left panels), 0.5 (center panels), and 1.0 (right panels). The 
upper panels show the case that the peculiar motion of sources is neglected, that is, the linear 
distortion and the nonlinear and the finger of god effects are omitted. The upper panels show 
the effect only from the geometric distortion. The middle panels show the case the geometric 
distortion and the linear distortion are incorporated but the nonlinear and the finger of god 
effects are omitted. The lower panels show the case all the distortion effects are incorporated. 
In this figure the ACDM model with the same cosmological parameters as those in figure 6 is 
assumed. 

Fig. 8 — Contour of ^2 C (-R)/^o C (-^) w hh fixed R = 10/i _1 Mpc for various cosmological models, 
(a) ACDM model with fi = 0.3, Ov = 0.7; (b) ACDM model with tt = 0.4, fi A = 0.6; (c) 
open CDM model with fi = 0.3, fi A = 0.0; (d) open CDM model with fi = 0.4, fi A = 0.0. 
In each model we adopted h = 0.7 and a = 0.5. The contour lines are —0.8, —0.5, —0.3 from 
left to right for the ACDM models. While the contour levels are —0.5, —0.3, —0.1 for the open 
CDM models from left to right. The region between the dashed lines in each panel satisfies 
4 < R c < 6. Here Ug is determined from the cluster abundance (Kitayama & Suto 1997). We 
assumed the cosmological redshift space with s = SEd(z). 

Fig. 9 — Same figure as figure 8 but with fixed R = 20/i~ 1 Mpc. The contour-levels drawn by 
solid lines are shown in each panel. 
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